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In low-scale gravity models, a particle collider with trans-Planckian collision energies can 

be an ideal place for producing black holes because a large amount of energy can be con- 

^N centrated at the collision point, which can ultimately lead to black hole formation. In this 

^^ article, the theoretical foundation for microscopic higher dimensional black holes is reviewed 

, and the possible production and detection at the LHC is described and critically examined. 
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I. INTRODUCTION 

A long standing problem in high energy physics is understanding the quantum nature of gravity. 
The problem is often defined in terms of the Planck scale, Mp, which characterizes the strength 
of the gravitational interaction and sets the cut-off scale of the effective quantum field theoretic 
formulations for the strong, weak nuclear and electromagnetic interactions. Over the past several 
decades, the Planck scale has been regarded as a fundamental scale of physics, but the highness 
of the scale has been also known to lead to the notorious hierarchy problem. The big hierarchy 
between the Planck scale (~ 10^^ GeV) and the electroweak scale (~ 10^ GeV) can be translated 
into the extreme ratio between the gravitational constant, Gjy = 1/Mp, and the Fermi coupling 
constant, Gp, which represents the relative weakness of the gravitational interaction compared to 
the weak nuclear interaction, as follows: 

On _ 6.709 x IQ-^^ GeV-^ _ ^^.3^ 



Gf 1.166 X 10^5 GeV-2 



r\j 



The presence of a high scale is problematic particularly in the scalar sector since there is a large 
quantum quadratic correction, Smf_i ~ A^. The Higgs sector of the standard model thus suffers 
from a fine-tuning problem at an extremely high level, 

. _ mLe + ^^Lioop (100 GeV)^ 34 

Atuning = J-, ~ (1019 GeV)2 ~ ^^ ' ^^' 

when A ~ Mp. It is evident that the level of fine tuning, Atuning 1 is essentially given by the ratio 
in Eq. M . The technical aspect of the fine-tuning problem is known to be resolved in the presence 
of additional symmetries such as supersymmetry and conformal symmetry provided that their 
breaking scales are close to the physical Higgs mass, which still demands additional (dynamical) 
explanations. 

In the past decade or so a new field has emerged in high energy physics and has expanded 
to other fields including quantum and classical gravity, cosmology, astrophysics and numerical 
relativity based on the realization that extra dimensions can lead to "low scale gravity" and "brane 
world" models. A higher dimensional brane world model with large, flat, extra dimensions (Arkani- 
Hamed, Dimopoulos and Dvah (ADD, 1998) [HE], Antoniadis+ADD (AADD, 1998) [3j) and a 
five-dimensional model with a part of highly warped AdSs space (RS (1999) |1]) have catalyzed an 
entirely new approach to the hierarchy problem. In these new higher dimensional models, a brane, 
like an ordinary domain wall, soliton or D-brane in string theory, is introduced to confine a subset or 



entire fields of the standard model. Gravitons, on the other hand, can propagate through the higher 
dimensional bulk and the effective strength of the gravitational interaction for a brane- localized field 
is diluted by the large volume of extra dimensions or red-shifted by the 'warp factor' along extra 
dimension(s). As a result, gravity becomes weak as has been observed in ordinary experiments. The 
true strength of gravity can only be measured at small distance scales below the compactification 
radius or the size of extra dimensions, but can cause interesting phenomenological consequences. 
On the true gravity scale. Mo, can be as low as the electroweak scale then the hierarchy problem 
does not exist, u 

The most striking prediction in low-scale gravity models is the production and decay of higher 
dimensional rotating black holes at the CERN Large Hadron Collider (LHC). Indeed, back in 1974 
Penrose first showed that a black hole can be formed in classical high energy head-on collisions 
[H] and in 1987 't Hooft [9] made preliminary arguments for graviton dominance and black hole 
formation in the trans-Planckian domain, which has been clarified by many subsequent studies. The 
production of a black hole during particle collision can be easily understood when we accept the hoop 
conjecture suggested by Thorne in 1972 |10j . The hoop conjecture essentially states that a black hole 
forms if and only if a large amount of energy is packed in a small region that can be surrounded by 
a hoop with the Schwarzschild radius of the energy. A high energy collider can provide an idealistic 
environment for this phenomenon as a large amount of energy can be concentrated at the moment 
of particle collision. In the models with large or highly warped extra dimensions, gravity becomes 
strong at a low scale in a TeV range so that one can conclude that a collider with a multi-TeV 
collision energy, such as the LHC, can produce black holes with TeV~ sized event horizons. Some 
rough estimation using the geometric cross-section for black hole production, a ~ 7rrs(-y/s)^, were 
made by Banks and Fischler (1999) [11], Giddings and Thomas (2001) |12] . and Dimopoulos and 
Landsberg (2001) |13| have drawn considerable attention not only from physics communities focused 
on high energy, string theory, general relativity, numerical relativity, and cosmology, but also from 
the general public including even opponents of the LHC experiment (see |14) and |15) for scientific 
responses to the opponents). We also note that Argyres et al. (1998) [16] considered the TeV scale 
primordial black holes as dark matter candidates or as seeds for early galaxy formation. 

After formation, a black hole radiates energy through Hawking radiation p!7kl9| mainly to the 
standard model particles on brane |2^. The decay signatures depend on greybody factors, s^lm.) 



^ To be fair, one should provide the mechanism to stabihze the 'radion' or geometry to completely solve the hierarchy 
problem. See Goldberger-Wise (1999) [S] for this mechanism. There is now a wide class of string scenarios, including 
the "flux compactifications" (2001) of Giddings, Kachru and Polchinski (2001) [6] and of Kachru, Kallosh, Linde, 
and Trivedi (2003) T and other general warped compactifications, which make the distinction between large volume 
and large warping less clear-cut. One can have either, or a combination of both, in more complete scenarios. 



and a thermal factor with angular velocity at the horizon Q, 

dEgPrn 1 s^em 



dtdoj 2-K e('^-"^^)/^ - (-1) 



2s' 



(3) 



where w, s, £, m denote the energy, spin and angular momentum quantum numbers of the radiated 
particle, respectively. In the semi-classical domain, greybody factors have been obtained for an 
arbitrary spin (s = 0, 1/2, 1) in |2Tti2i] . [25ti28] (also see |29H33] for various extensions). The full 
calculation for gravitational radiation {s = 2) is still not available, but there have been some recent 
developments |35H37| . Semi-realistic Monte-Carlo (MC) event generators for the black hole events 
at the LHC have been developed. The most developed MC generators are BlackMax [38J and 
CHARYBDIS |39l SO] (also see CATFISH [41j). BlackMax v.2 and CHARYBDIS v.2 are used for 
LHC black hole searches by the CMS and ATLAS collaborations. 

In this article, we review the theoretical and experimental foundations of this rapidly developing 
field. We begin in Section [IT] with an introduction to physics in higher dimensions, in which 
we explain low-scale gravity models with large and warped extra dimensions. We also give a 
short introduction to compact hyperbolic extra dimensions where light Kaluza-Klein states are not 



necessarily expected, but still the gravity scale can be low (~ TeV). Section HI describes the black 
hole solutions in higher dimensions where the black hole uniqueness theorem is not generally valid. 
Emphasis is placed on the Myers-Perry black hole. Black rings and multiple black hole solutions 



and their stability issues. Section IV is devoted to describing physics in the trans-Planckian energy 
domain. We explain how the known physics of combining quantum field theory and classical 
general relativity provides good descriptions in the trans-Planckian domain and review the black 
hole formation process by superimposing the two Aichelburg-Sexl solutions to estimate the cross 
section for making black holes in particle collisions. An extended phase diagram for the relation 
between the probing distance and the available energy is also explained at the end. In Section [V] 
we focus on the LHC experiment to identify black hole signatures. The expected signatures are 
obtained considering the precise greybody factors. Issues in the balding and Planck phases of the 
produced black hole are discussed. After reviewing the recent LHC result with ^/s = 7 TeV collision 
energy, we critically re-examine the interpretations. In Section [VT] we conclude with a discussion 
of future higher energy runs with y/s = 14(100) TeV and open problems. We present the relevant 
physical scales in Appendix [A] and briefly sketch the numerical method to solve the generalized 
Teukolsky equation to calculate the greybody factors in Appendix [Bj 

Some aspects of this topic have been described in other reviews by Landsberg |42j, Kanti 



Giddings |44| and Park |45| and also in lectures by Kanti [46]. 



II. PHYSICS IN HIGHER DIMENSIONS 

In this section we first provide a generalization of spacetime taking n-extra dimensions into 
account. Indeed, no particularly compelling theoretical reason exists for three spatial dimensions. 
One may note that only in (3 + 1) dimensions, the Kepler problem with the potential of the form 
V = —a/r allows closed orbits |48) and thus one may find an anthropic reason for the three existing 
large spatial dimensions. However, it still allows the existence of the extra microscopic dimensions, 
which are not responsible for the planetary motions. More importantly, consistency conditions for 
the known extension of Einstein gravity and the quantum field theoretical embedding of it, namely, 
string and supergravity theories, require a certain number of extra dimensions. Furthermore, it has 
been revealed that extra dimensions can provide totally new insights to reformulate the notorious 
old problems in particle physics including the big hierarchy problem, which will be discussed in this 
section. 



A. Einstein-Hilbert action in higher dimensions 

In (3 + 1) dimensions, pure gravity is described by Einstein-Hilbert action: 

SER = -ld'x^,^R, (4) 

where g = |Det(g(^i,)|. Setting the action dimensionless {h = 1, c = 1), the Newton's constant G has 
a dimension of [Length] = [Mass]" and thus defines a mass scale, l/yG. This scale is nothing 
but the Planck scale Mp = 1/VG = 1.22089(6) x 10^^ GeV [39]. 

In the presence of the extra n-dimensions, the Einstein-Hilbert action can be generalized to a 
form 

Seh = Si-^ Sd = - <fx VdD .r. ^ Rd, (5) 

where Gd has a dimension of [Length] ~ = [Mass]" as [Rd] = [Length]" . A higher dimen- 

sional generalization of Einstein equations could be obtained by varying the action in Eq. (|5|. 
One may define a generalized Planck scale in Z)-dimensions (see Appendix^ for Planck units in D 



dimensions): 



1 



where Nn is a numerical factor for which different definitions have been employed in the literature. 
Some representative definitions that are widely utilized include the Particle Data Group conven- 
tion (PDG) [49j, which is the same convention in Ida-Oda-Park (lOP) [21], the Randall-Sundrum 
convention (RS) [3], Giddings-Thomas convention (GT)[12j and Dimopoulos-Landsberg convention 

(DL) m-. 



Mn 



(27r)" PDG , lOP 

1 RS 

(7) 

2 • (2^)" GT 

Svr DL. 



In this review, we follow the PDG convention. Now, starting from the higher dimensional action, 
Sd, the four-dimensional action ^4 can be derived as an effective action after integrating out the 
extra coordinates. In this respect, Md should be regarded as a fundamental scale of gravity, from 
which the Planck scale in four dimensions is derived after compactification. 



Let us consider a pretty general model of an extra dimension with a "warp factor" P^(y*), which 
is a monotonic positive function of the extra coordinate y* . If we require the 4D Poincare invariance 
at each point of y*, the full metric of 4 -|- n-dimensional spacetime can be written in general as 

ds^ = W\y)j]^,dx^dx^ + 5n^n(y)d2/™dy^ (8) 

where we assume y* G [0,^]. We do not consider a negative warp factor value to avoid a naked 
singularity in extra dimensions. With the assumed monotonic property of the warp factor, we can 
link the change in the position along the extra dimension and the renormalization group "running" 
with respect to the energy scale. Thus one can impose the geometric meaning on the UV-scale and 
IR-scale in extra dimensions. 



B. Flat, factor izable geometry 



If the warp factor is trivial, W{y) = 1, and the extra dimension is flat, gmn = —^mm the resultant 
geometry is simple, thus we can easily compute the 4Z) effective action as follows: 



Sd 



'Jmn — '-'m 



W{yy- 



^- /"d" 



yV9y 



WttGd 



X / (fx^R + 



Vn 



IGttGd 



X j d'^x^R + 



(9) 



where we can find the relationship between the Newton's constant in 4D and in general dimensions: 



GrA 






(10) 



Here, the volume of the extra dimension is denoted hy Vn = j (F"yJ~g^. If the volume of the flat 
extra dimension is large, V^ — )• oo, we get a tiny Newton's constant G4 — )• 0, thus the huge Planck 
mass M4 —7- 00. The large extra dimensions provide a geometrical explanation for why we have 
weak gravity or a large Planck scale as was suggested by (A) ADD [H [3]. 

A particularly interesting possibility is that M/5 is as low as the electroweak scale of the order 
of 1 TeV (Af/) ~ ITeV) such that there is no hierarchy between the fundamental scale and the 
electroweak scale. The solution to the hierarchy problem works if 



Ml 



or 



Md^ (27r)" 
y„^/" ~ (27r) X 10^°/" X {1/Md) 



(MD'^Vn) ~ 10 
30 



30 



(11) 



10 n -Hm 



where we took M4 = I/VSttG ^ 2.4 x 10^^ TeV and 1/Md ~ 1/ TeV ~ 2 x lO"'' fm. The actual 
numbers for various dimensions are 



L = K^/" 



~ < 



10^2 m n = 1, 

1mm n = 2, 

10 nm n = 3, 

4 X 10^ fm n = 4, 

10^ fm n = 5, 

100 fm n = 6, 

20 fm n = 7. 



(12) 



Obviously the first case, n = 1, is ruled out by astronomical measurements. The second case, n = 2, 
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is also ruled out by the short-distance Torsion pendulum experiments carried out by Kapner et al. 
(2006) |50| and Adelberger et al. (2006) |51| . Due to the mass gap, (Am ~ l/T/^'") is small with a 
large volume, and Kaluza-Klein (KK) gravitons can contribute to astrophysical processes such that 
models with a larger number of extra dimensions (n > 3) can also be strongly constrained. The 
requirement that KK gravitons do not contribute to the supernova cooling process by carrying the 
energy away gives the bound Mo > a few TeV for n = 3 (Hanhart et al. 2000, 2001) |521 153] , An 
even stronger bound, M^ > 4 TeV was found by Hannestad-Raffelt (2001) [54j from measurements 
of diffuse 7 ray detection by the EGRET satellite for n = 3 and also for more involved processes 
with supernovas, neutron stars and pulsars (Hannestad-Raffelt (2002)) |55| . All these astrophysical 
bounds are for light KK gravitons Atti < 100 MeV, thus a nontrivial extra dimension model with 
a larger mass gap is not constrained by them. Taking the possible large astronomical uncertainties 
into account, n > 3 cases are still not excluded. 



C. Non-factorizable geometry 

When W{y) is a non-trivial function of y (again we still require a monotonic property), the 
geometry is warped and the effective theory has a scaling behavior, A(y) = W{y)A.. As in the flat 
case, we can read out the effective action by integrating the extra dimensions out as follows: 

Sd '-=-'-, - fd-yW\y)-^ X I d\^gW-^r''R,u + --- , 

W{y) J IQ-kGd J 

^" X / (fx^/^R + • • • , (13) 



IQ-kGd 

where barred quantities are canonically normalized. We can derive the conventional Newton con- 
stant out of the higher dimensional Gravitational constant as 

G4 = %^, (14) 

where "warped volume" is defined as 

Vn = j d''y^yW\y) (15) 

with the nontrivial warp factor W{y) in the metric along the extra dimensions. 
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For concreteness, we consider the n = 1 case in detail here 

m! = mI^, (16) 

which yields a natural relationships among dimensionful parameters, M4 ~ M5 ~ 1/V^i- For 
instance, an AdS space with W{y) = e~^ gives Vi = Jg dye~ ^ ~ 2k{^ ~ ^~ ) ~ 2fc '^here k is 
the AdS curvature. Thus, k ~ M5 ~ M4. 

Notably the standard model (or at least the Higgs field) is confined to the "IR-brane" located 
at y = i. (Conversely, the "UV-brane" is at y = 0.) The Lagrangian for the Higgs field is read as 
follows: 

^Higgs = / ^'^ / "^y^iy - ^W\v)C {W\y)r^^,, H) (17) 

= / (fxW\l)C {W\l)r,^,, H) . (18) 

By rescaling the kinetic term in C and taking the warp factor into account, the physical Higgs field 
with the canonically normalized kinetic term is defined as H-phys = W{l)H^ and the physical mass 
parameter on the IR (and UV) brane is given as 

Mphys(^) = W{ym = e)M, (19) 

Mphys(y/R) W{ym) ■ ^ ' 

In AdS space, Mphys = e~ M <C M such that the actual physical mass scale on the IR-brane is 
exponentially small compared to the one on the UV-brane. This phenomenon takes place due to 
the warp factor, which changes the "ruler" at each position y. Having exponential dependence, we 
can easily generate a big hierarchy between the UV-scale and the IR-scale as Randall and Sundrum 
(1999) jl] showed. A slice of AdS space arises as a solution to the five-dimensional Einstein equations 
with a negative cosmological constant (AdS curvature) in the bulk and well-balanced tensions for 
UV (IR)-branes. The setup provides an explanation as to why the IR-scale (or electroweak scale) 
is much lower than the UV scale (or Planck scale). In this setup, the fields on the IR-brane can 
"feel" the strong gravity on the IR scale Mir, ~ e~ M\j\r ~ e~ M5, and the Planck- TeV hierarchy 
is given by the warp factor: 

^^ ^ 10-15 = e-"'. (21) 

Muv ^ ^ 
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The most important lesson here is that gravity may become strong at the IR scale (~ TeV) 
in the presence of warped extra dimensions. If the extra dimension is large or highly warped, the 
strong gravity scale can be as low as the electroweak scale such that there is no hierarchy between 
the gravity scale and the electroweak scale. This is a unique feature of a higher dimensional setup 
that cannot be found in four-dimensional effective theories. This understanding opens up an entirely 
new area of extremely interesting phenomenological consequences including black hole production 
at the LHC. 

D. Compact hyperbolic space 

An intriguing possibility is that we live in higher dimensional spacetime, which contains compact 
hyperbolic extra dimensions. The compact hyperbolic spaces (CHS) have been relatively less inten- 
sively studied by physicists so far but its interesting mathematical properties have been explored 
since original work by distinguished mathematicians including Selberg(1965) |56) . Mostow(1968) 
|57) and Thurston (1978) ^58j. In particular, there are two properties listed below that deserve 
further attention. 

• Rigidity: Proved by Mostow (1968). 

• Large mass gap: Conjectured by Selberg (1965). 

A CHS can be obtained by taking a quotient space of a d— hyperbolic space {'H'^) with its funda- 
mental group (r): 'H'^/r. If we do not consider a cusped hyperbolic three-manifold or incomplete 
metric, the Mostow's rigidity theorem states that the geometry of a CHS of dimension greater than 
two {d > 2) is determined by its fundamental group. This means that once we fixed the local 
curvature of CHS, ic, and the volume, no more moduli remains to be fixed. This is largely different 
from other extra dimensions which often suffer from shape moduli stabilization issues. Of course, 
there remains the issue of super-selection of a particular fundamental group among many others. 

The large mass gap of CHS is based on the conjecture by Selberg (1965), who stated that there 
is a lower bound of the first eigenmode of the Laplace-Beltrami operator on CHS determined by 
the following inequality: 

kiic > \. (22) 

Indeed Luo and Rudnick (1995) [59] proved the conjecture up to kilc ^ ^171/784 ~ 0.22 using 
number theory approaches, and Cornish and Turok (1998) j60] and Cornish and Spergel (1999) |61| 
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also studied higher excitation modes on CHS. 

The presence of the lower bound opens up an entirely new perspective of extra dimensional 
physics. Usually, as in torus compactification, a large volume implies a small mass gap of the order 
of l/Vol„ '". Since the large volume is responsible for the small fundamental gravity scale (see Eq. 
(|12[)), a small mass gap or a light KK excitation of a bulk field seems unavoidable. However, due 
to the existence of the Selberg's bound, CHS does not have to accompany a light mode. This is 
a big advantage in terms of model building. If there is no light Kaluza-Klein state, the effect of 
extra dimensions is hidden and the model does not affect low energy experiments. In this case, the 
extra dimension is evident only through the strong gravity effect such as black hole production at 
particle colliders. This possibility was recently studied by Orlando and Park (2010) |62] (also see 
|63] for an inflationary model based on CHS). 

E. Summary and discussion 

In pure general relativity, the only input parameter is the number of spacetime dimensions. In 
string theory, a certain number of extra dimensions are required to be introduced to formulate the 
theory in a consistent way. Recently we have observed in many theoretical studies that the presence 
of extra dimensions provides new insights into many problems in particle physics too. 

On the other hand, we do not have any concrete experimental 'data' for extra dimensions thus 
we can only theoretically model the topological and geometrical structure of 'hidden' dimensions. 
Recently suggested 'large' and 'highly warped' and also 'hyperbolic' extra dimensions have opened a 
new window to addressing the problem of weak gravity or hierarchy in energy scales. A particularly 
interesting possibility is that the fundamental scale of gravity (which is supposed to be the largest 
scale of physics) may not be extremely large as was traditionally suggested by Planck energy, but 
can be as low as a TeV at least for a Higgs sector. If the low scale of gravity is really built in, we may 
learn about the presence of extra dimensions by experimentally probing the strong gravitational 
phenomena. The most dramatic phenomenon, which may not be that far from us, is microscopic 
black hole production and their decay at the LHC and sky ^\- Most probably those microscopic 
black holes are higher dimensional as long as their 'size' is much smaller than the compactification 
radius or local curvature radius of the background extra dimension: 

TBH < re. (23) 

However, one should note that the details of 'black hole signatures' depend on the detailed 



14 

structure of extra dimensions and also the precise distribution of matter fields in higher dimensions 
which we can only make models. Also we have to assume that our theoretical understanding of 
the quantum nature of gravity is limited. The best thing we can try now is making concrete and 
sound 'predictions' of a given set of models with theoretical confidence. A black hole can be well 
described by a semi-classical picture as long as its 'size' is big compared to the Planck length: 

io < TBH- (24) 

The next section discusses for such semi-classical and higher dimensional black holes. 



15 



III. BLACK OBJECTS IN HIGHER DIMENSIONS 

In this section, we review the black objects in higher dimensions. Indeed in higher dimensions, 
a generahzation of the black hole uniqueness theorem is not valid |64H69| . and a number of black 
objects can exist that share the same mass and angular momentum. For instance, the higher 
dimensional generalization of the Kerr black hole, or the Myers-Perry(MP) black hole, can share the 
same mass and angular momentum with black rings (BR) , the topology of which is distinguishable 
from the simple sphere of the MP black hole: S'^ x S^ (BR) vs. S"^ (MP bh) in five dimensions 
{D = 5). An even further set of solutions seem to exist in higher dimensions (D > 5), thus the 
classification of all stationary, asymptotically flat (or AdS, dS or hyperbolic) black hole solutions 
to the higher dimensional Einstein equation (in vacuum or with coupled matter) becomes a major 
problem in the study of higher dimensional gravity. 

In particle collisions, however, the most relevant solution is the MP black hole as other solutions 
seem to suffer from various instabilities (there is instability even for the MP black hole, especially 
at a high rotation limit). In this section, we focus on the MP black hole, black rings and other 
black hole solutions in higher dimensions. 



A. Schwarzschild-Tangherlini solution in higher dimensions 

The generalization of the Schwarzschild solution in higher dimensions is obtained by replacing 
1/r dependence of the Newtonian potential with 1/r as Tangherlini found in 1963 |70] : 

d4c^ = (l - ^) dt' - :p^^ - r'dnl_„ (25) 



r'- 



where a mass parameter fi is defined as 

where ri„ = i'Un+i)/2) ^^ ^^^ area of a unit n-sphere. Under the linearized gravitational perturba- 
tions, the D > 4 Schwarzschild solution is stable as Ishibashi and Kodama showed in 2003 |7HI72). 
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FIG. 1. The size of the Schwarzschild Black hole in D = 5, 6, • • • ,11 dimensions for a range of masses. One 
should notice that a black hole can be described by a classical GR only when r^ :^ l/Mjj such that one can 
neglect the quantum gravitational corrections. 



The event horizon locates at r '^ = fj,: 



rs = fJ- 



D-3 



IQ-kGdM \ D-3 



{D - 2)nD-2 



1 



ko- 



M \D-3 1 



Md 



Md' 



kn 



AfnTi^) \ "+^ 



+ 2)7r("+3)/2 



n 



(27) 



(28) 



where Mn is given in Eq. ([7]) for n(= D — A) extra dimensions. Explicitly, using the PDG convention 
where Mn = (27r)^, we found 



kn 



n = 0, 

1, 
2, 
3, 



(29) 



1/4tt ^ 0.08 

(2/37r)i/2 ^ 0.46 n 

(3/4)1/3 ~ 0.90 n 

(2/5^4) ~ 1.33 n 

(5^)1/5 ^ 1.73 n 

2(37r/7)^^^«2.1 n = 5, 

(1057rV2)i/'^ -2.4 n = 6, 

2(2^)1/73^/^ ~ 2.76 n = 7. 

Since the Schwarzschild radius increases as ~ M^'^^^^' the size slowly increases with respect to 
the mass in higher dimensions (see Fig. [I|. 
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The Hawking temperature and the horizon area (or entropy), respectively, are given as 

Th 
Sh 



Airrg 

1 dM r^-^nD-2 



D-2 dT 



4.Gd 



(30) 
(31) 



As M increases, r, and Sh also increase, but Th decreases. The smallest black hole is the hottest. 



'~h 
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FIG. 2. The thermodynamic variables (temperature, entropy) of the Schwarzschild black hole in D = 5 and 
D — 10 dimensions, respectively. 



B. Myers-Perry black holes 



In Z)-dimensions the number of independent rotation planes is A^ = [{D — l)/2]. In 4D, a unique 
rotation plane exists, but in 5D and 6D, for instance, there are two independent planes. In particle 
collisions, however, an angular momentum vector, which is perpendicular to the collision plane of 
the particles, is given so that a black hole solution with one angular momentum is the most relevant 
among all possible solutions. 

The higher dimensional rotating black hole solution was found by Myers and Perry (1986) |73] . 
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With one rotation, which is parametrized by a parameter a, the solution takes the form 

dslip = df - J^ {dt - as1d(t)f - —dr"^ - (r^ + a^)sld(l)'^ - r^c^df7|,_4, (32) 

where we used shortened notations sq = sin 0, cq = cos and 

^ = r'^ + a'^cl, (33) 

A = r^ + a^ - /i/r^-^ (34) 

The physical ADM mass (M) and angular momentum (J), which are obtained by comparing the 
metric at r = oo with the metric for a massive perturbation in Minkowski spacetime, are, respec- 
tively related to the /i and a parameters: 

J = ^^Ma, (36) 



where Go is Newton's constant in L)-dimensions (See Eq. ([7])). The a > parameter is the angular 
momentum per unit mass, which is common in the literature. The event horizon for D > 4 is set 
at A(r = th) = or 






and the area is given as 



AH = r\i-\ri + a')nD-2, (38) 



where a dimensionless rotation parameter a* = a/rn is recursively defined, r^ is the Schwarzschild 
radius for the corresponding mass M and the angular momentum J = 0. 



A higher dimensional Schwarzschild solution (i.e., spherically symmetric, static black hole solu- 
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tion) is recovered when a = 0. In this non-rotating hmit, 



E 


^r\ 






A 


^r^ 


- ^^|r^- 


-5 


J 


^0, 






f^ 


— r / ^ 





(39) 
(40) 
(41) 
(42) 



An upper bound on a exists for D = A where 



A(r;,) = rl- firh + a^ = {r^ - ^f + a^ - ^,^/A = 0, (43) 

a^ = /iV4 - (r^ - m/2)' < {^l/2f , (44) 

.-. a < /u/2. (45) 



Similarly for D = 5, 



A(rO = rl-fi + a^ = 0, (46) 

rl = f,-a^>0, (47) 

.'. a<^, (48) 



or jj3 < |l^ for D = 5. In contrast, for Z) > 5, there is no upper limit on a so that, in principle, 
a black hole can rotate with an arbitrarily large a. These highly rotating black holes are often 
called ultrasp inning, and are known to be unstable as shown by Emparan and Myers (2003) |74] . 
Recently, a bar mode instability was reported by Shibata and Yoshino (2010) |75] . 

The Hawking temperature and the angular velocity and the horizon area (or entropy), respec- 
tively, are given as 

^ {D-3) + {D-5)al D-2, 

_ 1 dM _ M 2al rf^HlD^ 

^^ - -w^w - W^Wh " i + ^J ^Gd ■ ^ ' 

where the limiting cases correspond to J —t- or a non-rotating limit. 
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C. Other higher dimensional black holes 

1. Black strings and black p-branes 

Since a direct product of Ricci fiat spaces is also Ricci flat, a new class of vacuum solutions can 
be constructed in D + p dimensions out of a vacuum black hole solutions with event horizons in D 
dimensions. Explicitly, a black p-brane solution is obtained by a direct product of the Schwarzschild 
black hole and a flat manifold: 

p 
ds%p-brane = ds^^h X JZ ^^dy'dy^ , (52) 

i 

in which the event horizon of the Schwarzschild black hole is extended in y* directions. Thus the 
horizon topology is from 7^ to ^ x R^. If the y* directions are compact by periodic condition y* ~ 
2/* + 27ri?j, the horizon topology becomes % x T^. lip = 1, in particular, the new black hole solution 
is called black string. Having an extended horizon, the black p-branes are not asymptotically flat. A 
'long' black string is known to be unstable as Gregory-Laflamme showed in 1993 |76| and 1994 |77| 
(also see Harmark, Niarchos and Obers (2007) [78] and Cardoso and Dias (2004) |79| for instability 
of a small Kerr-AdS black hole). 



2. Black rings 

In addition to the MP solutions, another black hole solution with horizon topology 5^ x S^ (or 
donut) was found by Emparan and Reall in 2002 [80] ■ The asymptotic geometry is flat. The S^ 
describes a circle that is stabilized by the centrifugal force of the rotating ring so that a lower limit 
on the rotating parameter exists different from the MP black hole. The solution is given in a form 

m 



F{x)' F{y) ^' {x-vY ' '\ F{y) ^ G{y) G{x) F{x) 

where F(x) = 1 + Ax, G{x) = {l-x'^){l + vx) and C = (A(A-i/)(l + A)/(l-A))i/2. The parameters 
A and u vary in the range 1 > A > z^ > 0, the coordinates —1 > y > — oo and 1 > x > —1. The 
axis of rotation around the ?/; direction is at y = —1. The other rotation axis around i;^ is x = 1 
and X = —1 inside and outside of the ring, respectively. The horizon is at y = — l/i/. 
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In five dimensions, the spin can also grow indefinitely, but only if the spinning object is a ring: 

where |^ is the upper limit of the MP black hole in D = 5. Due to the existence of the lower 
limit on the angular momentum, we may conclude that black rings may not be produced in particle 
collisions as discussed by Ida, Oda and Park (2003 )|21j. However the production of black rings is 
not conclusive in Z? > 5. 

3. Multi black holes 

In D = A there are no stationary vacuum multi-black hole solutions to the Einstein equations. 
However, such black holes exist in D = 5 and also in D > 5. These multi black hole configurations 
could not be obtained by initial data about colliding particles, but they have their own interesting 
issues (see however, Yoshino-Shiromizu |82) for black ring formation by closed string collisions). 

The simplest multi black hole solution in higher dimensions is 'black saturn', which was con- 
structed by Elvang and Figueras in 2007 [183]. Black saturn is composed of a Myers- Perry black 
hole in the center, which is geometrically surrounded by a black ring. The angular momentum of 
the ring makes two black holes remain apart and stay in balance. An interesting property of black 
saturn is that its total angular momentum measured at infinity can vanish. This can happen when 
the black ring and the central MP black hole rotate in opposite directions. This is very distinc- 
tive from the four-dimensional case for which the Schwarzschild solution is the unique solution in 
D = A with vanishing angular momentum (J = 0). Other more complicated configurations are also 
known such as the black di-ring discussed by Iguchi and Mishima (2007) [84j and bicycling black 
rings by Elvang and Rodriguez (2008) |85]. Multiple-MP black hole configurations have also been 
considered, but the existence of conical singularities seems to be unavoidable [86] . A good review 
of exact solutions for higher dimensional black holes was recently presented by Emparan and Reall 
(2008) |87| and Tomizawa and Ishihara (2011) 



D. Summary 

It is important to realize that the black hole uniqueness theorem is generally not valid. The 
different black objects in higher dimensions can share their masses and angular momenta even 
though their geometrical properties are completely different. However, among many new black 
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hole solutions that have been constructed in higher dimensions to date, the Myers-Perry(MP) 
and Schwartzschild-Tangherlini (ST) solutions seem to provide the best realistic descriptions for 
rotating and non-rotating black holes at the LHC. Indeed the MP and ST solutions have classical 
stability under linearized gravitational perturbations so that they are the most probable candidates 
for produced black objects in particle collisions. Keeping this fact in mind, we summarize the 
most relevant physical quantities for the ST black hole ( J = limit of MP black hole) in various 
dimensions for future reference. We introduce a convenient dimensionless parameter M = M/ M^, 
for the black hole mass M and explicitly write rs,Th and Sh in various dimensions: 



TsMd 



(2/37r)i/2(M/MD)^/2 « 0.46M1/2 „ ^ ^ 

(3/4)1/3(M/Md)^/^ n = 2 

(2/51/4)(m/Md)1/4 n = 3 

(57r)i/5(M/Mz))i/5 n = 4 

2(37r/7)^/^(M/MD)i/^ n = 5 

(1057rV2)i/^(M/Mz))i/7 n = 6 

_2(27r)V4/3i/S(M/MD)^/^ ~ 2.76M^/^ n = 7. 



(55) 



47rTh 
Md 



^/6^(M/Md)"^/2 ~ 4.34M-1/2 

62/3(M/Mz5)"^/3 

2 • 5^/^{M/Md)'^/^ 

{5^ / Try /^{M/Md)-^/^ 

(3^ • 7/^)i/6(M/Mz))-i/6 

(2 • 7^3 • 5 • 7r2)i/7(M/Mz))-i/^ 



1 
2 
3 

4 

n = 5 
n = 6 



n 



n 



n 



n 



^(26 . 3/7r2)i/8(M/MD)-^/^ ~ 2.90M-1/S n = 7. 



(56) 
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Y^3W27(M/Md)^/2 ~ 1.92M3/2 n = 1 

(3/4)2/37r(M/Mz))4/3 n = 2 

87r/5^/'^(M/MD)^/'' n = 3 

5/, = <( 2(57r6)i/5/3(M/Mz))^/5 n = 4 (57) 

8(3-7r77'^)i/6(M/Mz))'^/6 n = 5 

(105 • 7rV28)i/7(M/Mz))8A ^ = 6 

, • (27r5)V4/3i7/8(^/^^)9/8 ^ 3.85M9/8 n = 7. 
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IV. PHYSICS IN TRANS-PLANCKIAN ENERGY DOMAIN 

"The scattering process of two pointlike particles at CM energies in the order of Planck 
units or beyond is very well calculable using known laws of physics because graviton 
exchange dominates over all other interaction processes. At energies much higher than 
the Planck mass, black hole production sets in, accompanied by coherent emission of 
real gravitons. 
-G. 't Hooft (1987)" 

As 't Hooft provided these preliminary arguments in 1987 |9l [89], which have been clarified by 
many subsequent studies, gravity becomes dominant in the trans-Plan ckian domain and the 'known 
laws of physics', i.e., a combination of quantum field theory and general relativity provides a good 
description of certain features of trans-Planckian physics. In this section, we review the scattering 
problem in the trans-Planckian energy domain, which is described by the known physics in the 
't Hooft sense. We use the conventional terminology in particle physics ^/s instead of Ec.m. for 
the scattering energy in the center of the collision frame. Relatively little is known in the Planck 
domain (-^/s ~ M^), but a large part of the trans-Planckian domain {y/s ^ Md) and sub-Planckian 
domain (y^ <^ Md) can be understood in terms of known physics, i.e., perturbative quantum field 
theory and Einstein gravity in higher dimensions. 

A. The limits of the effective field theory approach 

In the effective field theory framework, an action for describing physics below a high scale or 
a cut-off scale A is composed of a set of quantum field operators, which may be or may not be 
suppressed by the A scale depending on their dimensions. An operator 04+^ of the mass dimension 
(4 -|- d) is suppressed by A as the Lagrangian density C has the mass dimension 4 or the number 
of spacetime dimensions. An operator that is suppressed by A'^'*'^ is regarded as unimportant or 
irrelevant at low energy, q'^ ^ A^, and its effect is regarded as a small correction to the effect not 
suppressed by A. The generic form of the Lagrangian density in flat spacetime with the Lorentzian 
metric Tj^i, = diag(l, —1, —1, —1) is given by the following schematic expression: 

£eff = O4 + O5/A + Oe/A' + • • • . (58) 

This approach clearly fails when one considers a physical process beyond the cut-off scale A. In this 
"trans-cut-off domain" all the higher order terms become equally important and their effects can 
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no longer be regarded as small corrections to the leading order result. A higher-scale new physics is 
needed before this effective field theory fails. A good historical example is Fermi's four fermi theory 
of weak interactions, which failed on high scales but the electroweak theory of SU{2) x U{1) takes 
over the role of the new gauge bosons W, Z and 7 (photon). 

Now we can consider an effective field theory for matter fields coupled with gravityjj A mass 
scale is set by the Planck scale, A ~ M4, so that the effective Lagrangian for matter fields Sm is 
given in terms of operators, which may or may not be suppressed by M4: 

S = Seh+ [ d^x^iOi + O5/M4 + Oe/Mi + •••)• (59) 

When one considers physical processes in a low energy regime with small energy transfers, q'^ <C -M4 , 
higher order terms 0>4 are irrelevant. On the other hand, when q^ ~ M|, all higher order terms 
become non-negligible thus this perturbative description of the effective theory would not work. 
Indeed, if the relevant energy is well below M4, gravity is weak (the weakest force of all known 
fundamental forces). One can consider the gravitational interaction a small perturbation to the 
known non-gravitational interactions. When the energy becomes high and reaches the Planck 
energy, the gravitational interaction becomes dominant over all other forces as the gravity directly 
couples to the energy-momentum tensor differently from gauge interactions. However, one should 
note that a large momentum transfer can be carried by many soft exchanged gravitons, in which 
case the higher-dimensional operators appear not to play a role [98]. 

In D = 4 + n dimensions, a higher order term in addition to the Einstein-Hilbert actions is 
written as 

S= f d^xy^^^- {R + O {R^/Ml) + •••). (60) 

The actual form of the second order term is suggested in string theory, e.g., aB? + hR^^'^R^^ -f- 
cR^^P" R^upa with {a,b,c) = (1,-4,1) or the Gauss-Bonnet term; there is no ghost in AD. The 
second and higher order terms can be neglected when the curvature (=Ricci scalar) is 'small' or 
R/M'^ <^ 1 where we get the conventional Einstein-Hilbert action. 



A variation of the Einstein-Hilbert action is given as SSeh ~ — ^g^g S d'^^y/O \R^i.v — \gfiuR\ Sg'^" . The energy- 
momentum tensor of the matter Lagrangian is defined by 5Sm = | / d'^^s/gTy.uSg^^ so that (JSeh + SSm = leads 
the Einstein equation Ry_^ — ^gf^„R — SnGTfj,^ — 0. 
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B. The trans-Planckian domain: ^s ^ Mo 

First let us consider relevant length scales for the trans-Planckian scattering of the two parti- 
cles. They are the Planck length {^d)i the de Broglie length of the colliding particles (^de)) the 
Schwarzschild-Tangherlini black hole length (r^) corresponding to the CM energy and the impact 
parameter of the scattering in terms of angular momentum |90j : 



fGDh\D-^ 
Airhc 



tdB 



1 _ 

(61) 

(62) 



-C.{^y-\ (63) 

24;, (64) 



1 
where Cd = ( (r)_o\Q ) ^^'^ ^D-2 is the area of the unit sphere, S^~'^. 



{D-2)Qd-2, 
It is interesting to notice in the limit with ?i — t- and Gd, ^/s fixed, 

Id -^ 0, ^dB ^ 0, r^ ^ finite, (65) 

which means that semi- classical physics is responsible for trans-Planckian scattering and the charac- 
teristic length scale for the dynamics is r^. For instance, when the impact parameter of the collision 
is 6, one can estimate the small scattering angle 9 ~ (rs/6)""'"^ as is confirmed by the Eikonal ap- 
proximation in the trans-Planckian regime (see earlier work in string theory |91H96) and also see the 
relatively recent work |90 | I97f[99] ) . As h approaches r^, non-linear gravitational attraction becomes 
important and at around & ~ r^, black hole formation is expected. 

A quantitative definition of the trans-Planckian domain [A) for a black hole is given in terms 
of the Kretschmann invariant /C = \R^^'^^ Rlmnp\^ / ({D - 2)^/{D - 1){D - 3)^|,) in D di- 
mensions by Okawa et al. using a convenient energy scale Ep = {{D — 2)ilj:)_2/16vrG£)) ~ 
[Tnni [TnT] as follows: 

trans-Planckian domain: A = {A : Inf /C > 1}, (66) 

where the inequality is saturated for a black hole with a horizon of the size of the corresponding 

Compton wave length of the same mass. In the regime A with Inf /C <^ 1, we can solely trust 

A' 
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classical approximation over which quantum effects are suppressed as 0{IC) <^ 1. 



1. Parton level collision 

When we consider the process of a two-particle collision in asymptotically flat spacetime we can 
define a reference frame for the center of mass (or CM frame) in which the total sum of spatial 
momenta vanishes at the asymptotic limit. The CM energy is also obtained by summing the energies 
in the same reference frame. If the CM energy is significantly larger than the Planck energy, we 
call the collision trans-Planckian or super-Planckian. If the particles participating in the collision 
are composite or non-elementary, the actual CM energy contributing to the collision should be 
obtained for 'partons', which could be quarks and gluons in hadron colliders such as the LHC. This 
CM energy of partons is called the parton level CM energy, which often is much greater than the 
sum of the rest mass energies of partons due to a large kinetic energy contribution of the partons. 
On the other hand, the parton level CM energy is smaller than the total collision energy because 
other particles (spectators) do not contribute to the collision, but still carry energy away. If two 
partons have the same mass, the parton level CM energy is 

EcM = ^^^^== = 2m-f{v), (67) 

where t; < 1 is the velocity of each particle in the CM frame. As defined above, the collision is 
called trans-Planckian when -EcM ^ M^ in D-dimensions or 

2m^{v) ^ Md trans-Planckian, (68) 

thus for partons of mass m <^ M^, a large 7 factor is required. For instance, light quarks (up-quark, 
down-quark) and gluons are partons in a proton, and each parton carries a part of the proton energy 
(E'parton = a^-E'proton ~ 7(a^?72pi.oton)) < X < 1) but still Can be larger than its rest mass energy. The 
required 7 factor is estimated as 

7(^)»^=^^^^=~^:^^^^~^^(MB/TeV)xl03, (69) 

'XiX2'ml^^^^^ V^ia;2 "Iproton V^ia:2 

which defines the ballpark parameter range for the LHC study. On the other hand, when we consider 
'heavy' objects, especially colliding black holes, their masses could be assumed to be much larger 
than Md so that we may still reach the trans-Planckian domain with a rather small 7 ^ C(l)- 
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This domain is nothing but the classical physics domain. 



2. Black hole formation 




FIG. 3. Scattering trajectories of a slowly moving particle by a stationary particle, a slowly moving particle 
and an ultra relativistic particle, respectively (from left to right). The black thick lines are the trajectories. 
The colorful thin lines indicate the gravitational field associated with the target particle a.t u — z — vt — 0. 



A particularly interesting phenomenon is black hole formation in trans-Planckian collision with 
a sufficiently small impact parameter. The black hole formation in particle collision can be easily 
understood by considering two infinitely boosted particles (take m — )• 0, w — )• 1 and y^ = fixed). 
When a particle is highly boosted, the field line of gravity of the boosted particle is well confined 
on the transverse plane of the momentum. This highly boosted gravitational field is described by a 
gravitational shock wave identified by Aichelburg and Sexl |102) and the superimposed solution by 
Dray and 't Hooft in 1985 [89] for £> = 4 and Ferrai et al. for L> > 4 in 1988 [TOS]. (see [M] for an 
intuitive exposition on black hole formation in high energy collisions). In Fig. p] for illustration, 
we plotted the gravitational field lines for a rest particle, a boosted particle and a highly boosted 
particle, respectively, from the left to the right. The solid black curves describe the orbit of a slowly 
moving particle in the presence of the boosted gravitational field. A slowly moving particle can 
'feel' the gravitational attraction at a distance if the source particle is at rest (figure on the left). 
However, if the source is highly boosted (the right-most figure), because the field is confined on the 
transverse plane, the slowly moving particle feels the gravitational attraction right at the moment 
of the encounter with the gravitational shock wave not before and after. Thus the slowly moving 
particle behaves just like a free particle before and after the encounter (see the third figure in Fig. 
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^. On the transverse plane, the gravitational field satisfies the Poisson equation 

V2$(x^) = -WTrGof^S^'^ix'), (70) 

where x* denotes transverse coordinates. The solution is easily found: 

( 16TrGDfi/{nD-3{D - 4)(x^Xi)(^-2)/2) ^D = A + n. 
The metric is now given as 

ds^ = dudv - dx'^ - ^{x')6{u)du^, (72) 

where the light cone coordinates are u = t — z, v = t + z and r and (p are the coordinates on the 
{D — 2)-plane. The metric is flat except for u = Orl 

In Fig. [4j we plotted a spacetime diagram of colliding shock waves. A shock wave traveling 
along v = t + z = Om the —z direction will not be influenced until the shock waves collide so that 
the superposition of the two shock wave solutions gives the exact geometry in regions I, II and III 
outside the future light cone (region IV) of the collision of the shocks. 

Having the explicit combined shocks, one can check the appearance of the apparent horizon 
(AH, a marginally trapped surface S, the outer null normals of which have zero convergence) in 
regions I, II and III for various values of impact parameter b. The best bound on the maximum 
impact parameter, 6maxi can be obtained on the Yoshino-Rychkov slice (outside region IV). If the 
apparent horizon (AH) exists, it is necessarily inside of the event horizon so the presence of AH 
guarantees the formation of an event horizon (EH) or a black hole according to the area theorem. 
Penrose (D = 4, 6 = 0) [H], Eardley-Giddings (Z? > 4, 6 > 0, a trapped-surface construction) |105] . 
Yoshino-Nambu (Z) > 4, 6 > 0, numerical solutions) |106) and Yoshino-Rychkov (Z? > 4, 6 > 0, 
refined numerical solutions) |107| studied the formation of AH in shock collisions for various cases. 
From 6max one can directly calculate the lower bound of the cross-section of black hole formation 



^max/ 



3 



Adopting a convenient unit of length ro = (STrGo/i/nD-a)^ and the continuous coordinates, 

u —^ u, (73) 

V - 2logre{u) + ue{u) ,D = 4: 

V + 2e(u)/(D - 4)r^-* + ue{u)r^-^° ,D = 4 + n. 



(74) 



r I 1 — itl 



'(w)^'"°) , (75) 

(P^^(P^, (76) 

the metric is given 

ds"" = dudv - (l + (D - 3)^^e{u)y dr"" + r' (l - -^e{u)y dQl_3, (77) 
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t-0 



z + t^O 



FIG. 4. Spacetime diagram denoting two shock wave collisions. The left parton moves along z — t — with 
the speed of light and right partons move in opposite direction and collide at the center. Regions I. II and 
III are described by the combined Aichelburg-Sexl (AS) metric and are thus completely determined. Right 
outside of the Yoshino-Rychkov (YR) slice (outside of IV) the AH formation is numerically determined, and 
thus the bound for 6Max can be calculated. The future light cone (IV) is not determined by the AS metric 
such that ^Max on the YR slice is the best lower bound on ^Max- 



The inaximuni value of the impact parameter can also be obtained under the assumptions of 
the Hoop conjecture, which was originally suggested by Thorne in 1972 [lOj. In four dimension, 
the Hoop conjecture states that a black hole with horizons forms when and only when a mass 
M is compacted into a region with a circumference in every direction, C < AttGaM, where G4 is 
the Newton's constant in four-dimensional general relativity. For higher dimensional black holes, 
Ida and Nakao [108j suggested an isoperimetric inequality Vd-3 ^ G^M, where Vd-3 denotes 
the volume of the typical closed (D — 3)-section of the horizon (hyperhoop), which is apparently 
different from a naive expectation C < (GM)^''^"^'. 



For colliding particles, Park(2001) |110j and Ida-Oda-Park(2002) [2lj generalized the original 
Hoop conjecture by taking into account the angular momentum as follows (also see work by Yoo et 
al.(2009) |109| ). Let us consider a collision of two massless particles with a finite impact parameter 
b and CM energy -y/s such that each particle has an initial energy y/sjl in the CM frame. The 
initial angular momentum is also determined by Ji = 6-y/s/2 in the CM frame. Suppose that a 
black hole forms when the two particles come into a region surrounded by a Hoop with a radius of 
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FIG. 5. Schematic view of black hole formation in particle colUsions. The Hoop conjecture provides a 
criterion for a black hole to form: b < 2rh{M, J). Here ^Max saturating the inequality should be understood 
as the infimum of the maximum value for the allowed impact parameter for black hole formation. 



the event horizon of the mass y^ and angular momentum J. Then the infimum of the maximum 
impact parameter for a black hole to be formed during the collision is 



Inf [Max(6(Vi)] = 2rhiVs, J = h^sjl). 



(78) 



Hereafter, for the brevity of notation, we simply denote 6Max = Inf [Max(6(-ys)]. 



Using the relations in Eq. (36) and also in Eq. (37), one can easily found the maximum value 



of the impact parameter is set saturating the the inequality h < buax [HI I110| : 

2 



tMax(\/s) 



1 '^ 
D-2^2^15-3 



iV~s), 



(79) 



(l+(^)') 



which can be compared with the numerically obtained values by Yoshino and Rychkov. When b = 
&Max) the rotation parameter a* takes the maximal value (a*)Max = {n + 2)/2 and the corresponding 
angular momentum is JMax = ^Max-v/s/^- Assuming the production rate is unity when b < ^Max; 
the production cross section of black hole can be expressed as 



cr{Vs) > vr6Max(\/s)- 



(80) 
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In Table IV B 2 we summarize the results obtained by Yosliino-Rychkov and Ida-Oda-Park and 
compare them with the naive value obtained by dimensional analysis. One should notice that the 
analytic calculation in Ida-Oda-Park assumes (1) no attraction of the geodesies before the collision, 
(2) no radiation of energy during the collision; both of these assumptions reduce the estimated 
maximum value of b or enlarge the ratio a/irb'^g^^. In this sense, the results by reported by Yosliino- 
Rychkov are remarkably consistent with the results described by Ida-Oda-Park. 

It is possible to study the perturbation of the geometry by superimposing two equal Aichelburg- 
Sexl shock waves traveling in opposite directions in D = 4 + n dimensions. D'Eath and Payne 
(1992) studied black hole-black hole collisions in D = 4 and made a crude estimate of the efficiency 
of gravitational wave generation, about 16.4% |lllflll3] . More recently, it was found that about 
25.0% {D = 4) and 41.1% {D = 10) of the initial energy is radiated as gravitational waves in 
head-on collisions |114| . It is interesting to note that the percentage of the emission increases 
monotonically with respect to the larger number of dimensions. This result is understandable as 
the number of degrees of freedom in gravitational waves in higher dimensions, is bigger than the 
one in lower dimensions. 



D 



o'yr/o'iop 
aYR/Trrs{Mf 



6 7 8 9 10 11 



0.71 1.25 1.57 1.69 1.73 1.75 1.73 1.72 
0.71 1.54 2.15 2.52 2.77 2.95 3.09 3.20 



TABLE I. The value of ctym obtained by Yoshino-Rychkov [107J in comparison with the value analytically 
obtained by Hoop conjecture, ctiop = ''''^maxi ^^ Ida-Oda-Park |21] and also a naive dimensional analysis 



max' 



3. Numerical relativity on high velocity collisions 

Only a limited amount of information could be obtained by Aichelburg-Sexl solutions. One may 
guess that a black hole forms when b < 6max such that a differential cross-section of black hole 
formation is simply given by a geometric cross section: 



da 
lb 



Y' (»^' 



or da/dJ oc J since J cc b [21], which gives the conventional result a = Jn^''^ db-^ = '^b'^^^. 
However, the actual formation is dynamical and the process is violent. Numerical relativity is one 
of the tools used to study this dynamical situation. 
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For the LHC study, the following quantities are absolutely important for a given set of initial 
data {y/s, b): 

• How much energy will be lost during the high energy collision? (Graviational radiation during 
the balding phase): -Bfinal 

• How large of an angular momentum will be lost during the collision? (Angular momentum 
extraction during the balding phase): Jfinal 

• Production cross-section for a given {E, J)finai: dEdJ 

Indeed we have observed impressive progress in this direction during the last five years or so. 



The ballpark value for j{v) > 10 (see Eq. (69)), which may not be reachable within a short 



foreseeable time scale but still the recent progress could provide valuable information. 

• Higher dimensional numerical relativity formalism: Sorkin (2009) |115) . Shibata-Yoshino 
(2009) [Il6], Zilhao et al. (2010) [HZl, Lehner-Pretorius (2010) [TT8] . 

• Numerical simulation for high velocity collision: Yoshino et al. |119) . Sperhake et al. (2008) 
[EO], Shibata et al. (2008) [m], Sperhake et al. (2009) fm\, Witek et al. (2010) [ 123] . 
Okawaet al. (2011) pHO] . 

• Bar-mode instability: Shibata and Yoshino (2010) |124] . 

4-. Planck domain: ^/s ^ Md 

Relatively less is known for this domain. In this domain, the de Broglie wave length is comparable 
to the Planck length so that a geometrical notion of length does not have a concrete physical sense: 

Ade Broglie ~ -^Planck- (82) 

In string theory, if the string scale is below the Planck scale, Ms ~ TeV < M^ "stringy effects" could 
show up in the Planck domain. In principle, string-black hole correspondence, stringy excitation 
modes of the order of Ms, string balls as a highly excited string state, and the a characteristic 
high energy behavior of string scattering cross-sections could be included in the list of "stringy 
effects" and some quantitative predictions have been attempted (see Cullen, Perelstein, Peskin 
(2000) Horowitz(1996) |125], [l26], Dimopoulos (2001) [127], Oda-Okada (2001) [l28], Anchordoqui 
et al.(2009) [l29], Nayak(2009) [T30] V 



C. A phase diagram 
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FIG. 6. A phase diagram describing the minimum distance scale to be probed with a given energy (momen- 
tum transfer) . Below the Planck energy AId , the minimum distance is determined by the uncertainty relation 
Ax > h/Ap. Above a certain energy, presumably Planck energy, the information behind the event horizon is 
hidden and the minimum distance is essentially determined by the size of the classical Schwarzschild radius 
''Sch ^ i<;i/("+i) for the given energy. The domain with a large distance Ax '3> h/Ap corresponds to a 
classical domain where quantum corrections are small. Ax < h/Ap is unphysical. 



A suggested 'phase diagram' in Fig. p\ explains the energy-distance relation (see Park (2008) 
, Giddings (2010) |98|, Dvali et al. (2011) [131J ). In the figure, the horizontal axis is for 
the amount of available energy in physical processes such as the scattering of two particles. The 
vertical axis represents the smallest distance scale that can be probed by the corresponding energy. 
In a domain well below the Planck energy (M^i), the distance scale is limited by the uncertainty 
principle or Ax > 1/Ap. As the energy increases, the quantum gravitational effects become more 
and more relevant and finally dominate over all other interactions. On the other hand, once the 
energy surpasses the Planck domain, a new window of black hole production would be available and 
then the smaller scale physics {i <^ Id) would be hidden behind the event horizon corresponding 
to the energy such that finally the smallest distance scale is essentially given by the size of the 
Schwarzschild radius r^ ~ Gs. In the presence of extra dimensions, the energy dependence of the 
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Schwarzschild radius is r^ oc M^'^^~^' in D-dimensions (more discussion is in the next section) so 
that Ax > Ap^/(^"-'^). 
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V. BLACK HOLES AT THE LHC 

The primary goal of the LHC is to find the Higgs boson to confirm the standard model of particle 
physics. The other goal, which is equally important, is to search for signatures from new physics 
around 1 TeV. Low scale gravity models are subject to be tested by the LHC. The various quantum 
gravitational effects in particle interactions as well as the black hole formation are evident in these 
models. Currently, the LHC accumulates data with a high integrated luminosity and ^/s = 7 TeV 
collision energy in pp-hame. Before the scheduled shut-down at the end of 2012, the accumulated 
luminosity is expected to be more than lOfb^ . At full power, the LHC is designed to produce 
luminosity greater than 10 — lOOfb" per year with a center of collision energy of 14 TeV. 

For the black hole, the basic observable quantity is the number of specific "black hole signals", 
which can be expressed as: 

A^BH signal = J^LHc X a{pp ^ BR + X) X 5r(BH ^ signal), (83) 

where £lhc is the luminosity of pp collision, which is accumulated by the LHC at a given time, 
a{pp —7- BH + X) is the production cross-section for BH production and i?r(BH — )• signal) is the 
branching ratio of the decaying black hole to the signal. 

The production cross-section of BH depends on the collision energy as previously discussed. The 
greater the collision energy is, the larger the cross-section: cr{E) ~ (GdE)'^'^ ' for semi-classical 
black holes. Thus, a higher energy run with a larger collision energy (-^/s = 14 TeV) will have an 
obviously better chance of producing black holes than the current 7 TeV run. If possible, a future 
high energy collider with even higher CM energy, e.g., the VLHC with 100 TeV CM, is definitely 
more desirable. Giddings and Thomas (2001) |12j and Dimopouplos and Landsberg (2001) |13j 
estimated A^bh and concluded that the black hole production rate could be large, with Md ~ 1 
TeV, but later developments also showed that there could be significant suppression of the rate. 
The most important factor is the threshold of black hole production. For a thermal or classical 
black hole, the energy threshold should be much larger than M^ or y/s ^ Mo, otherwise quantum 
gravitational effects easily spoil the theoretical predictions. Based on the entropy counting, for 
instance, the threshold should be set to at least Mmin/M/) > 5 for semi-classical black holes. With 
this threshold, the cross-section is largely suppressed by PDF as shown by Cheung |160| . 

After a black hole is formed, it decays and emits Hawking radiation. The actual decay products 
depend on the number of extra dimensions, the mass, charge and angular momentum of the black 
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hole, the wave function profile (the location) of particles in extra dimensions and also the geometry 
of the extra dimension. However, the large entropy of a black hole and the basic property of 
Hawking radiation determine some robust properties of the signature: 

• Large multiplicity 

• Flavor blindness 

The large multiplicity of the signal is inherited from the large entropy of the thermal black 
hole. Thermal black holes typically have a large mass (^ M^) and the Hawking radiation contains 
large numbers of particles, jets, in particular. The flavor blindness of the black hole signal can be 
understood since the Hawking radiation is essentially thermal. Even when we take the greybody 
factors into account (which we will discuss later), flavor blindness remains. Statistically we would 
expect the same number of electron, muon and tau particles in the Hawking radiation. However, 
Hawking radiation depends on the spin of the particle and also on the number of degrees of freedom 
for each particle. For instance, a gluon has eight colorful degrees and the spin is one. An electron, 
on the other hand, does not carry any color charge and its spin is 1/2 such that the portion in 
Hawking radiation is relatively smaller than a gluon. 



A. The production cross section in a pp collision 

The LHC is a proton-proton {pp) collision machine. Each proton is composed of several partons, 
and each parton carries a part of the proton energy and interacts with other partons or passes by as 
a spectator at the moment of collision. The actual collision energy therefore is less than the total 
energy of protons, but a part of it. Fig. ([7]) is a schematic of pp collision that produces a black hole 
of mass ^/xlX2S and angular momentum J. Considering the parton distribution function (PDF), 
and the fi{x, Q^) of each parton (i) |132| 1133) . the hadron level cross-section is obtained as: 

app^BU+x = du — V'/i(u,s)/j(-,s)<Tij^BH, (84) 

where o"jj_>bh is the cross-section of the two partons (z and j). The sum runs i,j = {q,q,g) and 
q = (u, d, s, c, b), i.e., partons in a proton. 
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FIG. 7. A schematic view of black hole formation in pp collision. 



1. Form factor 



The parton level cross-section of black hole production is formally given as a geometric cross- 
section modified by a "form factor" J^: 



^ab^BU = F-KrsiVsf 



(85) 



where r^ is the Schwarzschild radius of the given CM energy ys (see Eqs. (28) and (29) for explicit 



expressions). Here the hatted variables are parton level quantities. In general, ys, the collision 
energy provides the upper limit of the produced black hole: M < \/s considering the possible energy 
loss before and during the formation process. The form factor J^ implicitly includes information 
about energy loss, angular momentum, charge and other dynamical modifications. 

A commonly used expression for the form factor is obtained following the relation given in Eq. 



(79): 



J^lOP 



1 + 



' D-2-' 



D-3 



(86) 



This form factor describes the geometric cross-section of a 'black disk' taking the angular momentum 
of the black hole into account |2H 1110) . 

Another commonly used form factor is from a numerical work by Yoshino-Rychkov (2005) |107| 
using the Aichelburg-Sexl solutions discussed in the previous section. Here we provide a useful 
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FIG. 8. The fitted function for YR and lOP. The dotted points are numerical values in YR and the solid 



curve is the fitting curve |J^ = (1.75451 - 27.7194e"°-8i™'^^) in 4 < Z? < 11. 



approximate formula to fit J^yr to J-iop (see Fig. [8| : 



-^YR= {ai-a2e "^^)-Fiop, 
(ai,a2,a3)fit = (1.75451,27.7194,0.817237). 



(87) 



Assuming Vs = 10 TeV, Md = 1 TeV and J^ = J-jop as reference inputs, the parton level 
cross-section for D = 5 — 10 dimensions is roughly a ss (3.2, 6.4, 10.3, 14.7, 19.7, 25.1) nb. 



2. Split fermion, spin, charge and other factors 



One should note that the parton level cross-section can be largely modified if fermions are split 
in extra dimensions |134) . The splitting is a unique phenomenon in extra dimension models and it 
may be responsible for the observed Yukawa hierarchy |135| (also see |136l 1137) for universal extra 
dimension (UED) setup) and the proton longevity problem (Adams (2000) [139j) in particle physics. 
In Fig. |9] we drew three fermions in an extra dimension with splitting. The wave function of /i has 
a relatively larger overlap with /2 than with f-^. The location of the wave function is controlled by 
the bulk mass. The wave functions /i and /2 are relatively closer, but /a is a bit more separated. 
As a result, the wave function overlap between /i and /2 is larger than that between /i^2 and /s- 

To take the splitting into account, Dai et al (2006) |138] introduced an extra dimensional weight- 
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FIG. 9. Split fermions in extra dimensions, /i and /2 are rather closed located but fi, is at a distance. With 
this splitting, the probability of making a black hole is weighted by the wave function overlap integral. 



ing parameter w^^{b) for a given impact parameter b on 'our brane': 



I V max " 

W'\b)^ / P''{be)dbe 

JO 



(88) 



where P^^{be) is the probability that the particles' separation in the extra dimensions will be be- 
With the weighting parameter, the total cross-section can be expressed as 



^JMnVs Ju V Jo 



1 /"^ A; /■''Max 

du — dbw'^27:bfi{v,Q)fJu/v,Q), 

Md^s Ju V Jo 



(89) 



which can be significantly suppressed by small overlaps among split fermions. The splitting can 
also be important in warped extra dimension (see e.g. Rizzo (2006) |140| ). 



Hoop conjecture does not indicate whether the spin, charge and other quantum numbers of a 
collapsing or colliding particle may cause differences in black hole formation, but we at least know 
that on some dynamical time scale they may cause differences. Yoshino and Mann (2006) |141| and 
Gingrich (2007) [142J studied the effect of charged partons on black hole production and concluded 
that the black hole cross-section can be significantly reduced (also see Yoshino et.al. (2007) |143j ). 
In most papers on black holes at the LHC, it is assumed that the extra dimension is asymptotically 
flat and large so that the Schwarzschild-Tangherlini or Myers-Perry black hole provide a reasonably 
good description of the black hole. However, 'brane' may carry non-vanishing tension and the extra 
dimension may have a finite curvature in the vicinity of the produced black hole. These effects can 
also affect the production rate, but quantitative details are not known. 
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B. Black hole decay by Hawking radiation 

In the first stage of production, the concentrated energy is not settled down to a simple, sym- 
metric black hole, but rather violently emitting gravitational waves. This stage is called the balding 
phase after the 'no-hair' theorem. The spin-down phase follows after the balding phase. In the 
spin-down phase, the black hole is well described by the MP solution with one rotation parameter. 
In this stage. Hawking radiation takes energy and angular momentum away so that the rotation 
of the black hole slows down and proceeds to the non-rotating phase or the Schwarzschild phase. 
In the Schwarzschild phase, the rest of the energy is emitted as Hawking radiation, which is rather 
symmetric and democratic (equal for each species). At the end of the Schwarzschild phase, the 
mass of the black hole becomes not much greater than Mq for which a full quantum gravitational 
description is required. We schematically illustrate the evolution of the black hole produced in 



particle collision in Fig. V B| including the balding, spin-down, Schwarzschild and the final Planck 



phase where M ~ Mj:,. The phases were first described in [12j . 

1. Black hole evolution: Balding, spin-down and the Schwarzschild phase 

A black hole, that is first formed in a particle collision is expected to be highly asymmetrical. 
Considering the violent environment of colliding and merging particles, no particular symmetry is 
expected for a black hole in the early stage of production. Thus, in the balding phase, the black 
hole loses its "hair" by emitting energy, angular momentum and charge in the form of gravitational 
radiation and gauge field emission. At the end of this phase, it will settle down to a stationary, 
rotating black hole then it can be nicely described by the Myers-Perry solution. An upper bound 
for the angular momentum was estimated by taking the initial condition of the angular momentum 
into account using the hoop conjecture |110j : < a = (n -|- 2) J / {2M^hR\)h) < {n-\- 2)/2, but there 
is a much room for studies aiming toward understanding the balding phase. One may hope that 
numerical studies for dynamical black hole formation can achieve an improved understanding in 



the future (see Sec. IV B 3). 



The phase following the balding phase is the spin-down phase. During the spin-down phase, 
the black hole becomes more symmetric as it loses its angular momentum. As seen in Fig. [TT] the 
apparent horizon geometry becomes more spherically symmetric by emitting angular momentum 
and proceeds to the Schwarzschild phase. 

In spin-down and the Schwarzschild phases, the black hole loses its energy and angular momen- 
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FIG. 10. A black hole evolves via four steps-balding, spin-down, and the Schwarzschild and Planck phases. 
The balding phase is a dynamical stage of black hole formation during which gravitational radiation takes 
the energy and angular momentum away, but the precise rates are still unknown. The spin-down phase and 
Schwarzschild phase are essentially described by Hawking radiation, then the final Planck phase meets the 
quantum gravitational uncertainty. 



turn (when J > 0) through Hawking radiation |18] : 



dtduj V J 



where 






jm 



e^/n - (-1) 



2s 



U) 



m 



(90) 



OJ 



n 



UJ — TTlQ,, 

(l + a2)r/.' 



(91) 
(92) 
(93) 



are the reduced energy and the angular velocity on the horizon, respectively. The numerical factor 
sF^ni is called the greyhody factor (also known as the transmission coefficient), which describes the 
modification of the radiation in the vicinity of the horizon due to the curved geometry. Indeed, 
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FIG. 11. Evolution of black hole via Hawking radiation. The geometry becomes more spherically symmetric 
in spin-down phase. When the angular momentum is lost completely, Schwarzschild phase takes over then 
black hole evolves toward the Planck phase. 



taking the greybody factor into account, the spectrum of Hawking radiation to be measured by an 
observer at infinity does not look quite black-body like. The greybody factor for a given spin-s of 
the radiating particle is defined to be the rate of absorption by the hole for a steady in-falling flux 
at infinity. 



AT™ — Ar°"* 



.N; 



(94) 



im 



under the boundary condition that the radial wave is purely ingoing at the horizon. The greybody 



factor can be found by solving the generalized Teukolsky equation in Eq. ( 95 ) and the spin- weighted 



spheroidal harmonics in Eq. (96). 



The generalized Teukolsky equation in D > A was first derived by Ida-Oda-Park (2002) 
utilizing the Newman- Penrose formalism developed in |146| I147| . The equation turned out to be 
separable as the projected MP metric is Petrov type D. For a D = 4 Kerr black hole, the original 
Teukolsky equation was derived and solved in the early 1970s by Teukolsky 1 148m 50] and Page 
|15H 1152) . A recent work demonstrated that exact solutions for the D = 4 Teukolsky equation can 
be found (see, e.g., Fiziev (2010) [153j ). The radial equation of the generalized Teukolsky equation 
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FIG. 12. Greybody factors for a fermion field in D = 5. 



for D = A + n dimension |21j is 

dr \ dr J 

K"^ / K(2r + (n-l)nr-'^) 

+ -—+s{4iujr-i — ^ ^ — 

A V A 



(95) 



n{n - l)fir " + 2mau} - a lo - A 



R = 0, 



where A(r) = r'^ + a^ — fir'"^^^, K = (r^ + a?')uj — ma and A is the angular eigenvalue, which 
is an integration constant coming from the angular equation. The angular equation is called the 
spin-weighted spheroidal harmonic equation, which is common for D = A and Z? = 4 + n (see 
Goldberg et al. (1966)[E1], Leahy, Unruh(1979) [IM], Leaver(1985)[156j, Seidel (1988) [EI] and 
Casals et al. (2009) [T58]): 



1 d f . dS 



sinOde V de 



smi 



j + [(s - aw cos ef -{scote + m esc Of - s{s - I) + A] S = 0. (96) 



At low energy expansion, in the regime corh <C 1, the analytic form of the greybody factors for 
fields with spin-s = 0, 1/2, 1 were found by Ida-Oda-Park (2002) |21| using the 'matching method': 
the radial equation is solved near the horizon and at infinity then matched in the overlapping 
region. The matching condition determines the wave function in the whole regime in r/r/^ G (1, oo). 
(For a non- rotating black hole in D > 4 at the low energy limit of the radiation see |144| I145| ). 
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FIG. 13. Greybody factors for D = 5, vector field. Super-radiance modes (i.e., F < 0) are indicated. 



However at high energy, urh > 1 , the matching method does not work so that one has to rely on 
numerical methods. The first numerical solution for a scalar (s = 0) for the whole energy regime, 
cjrfi £ (0,oo), was found by Ida-Oda-Park (2004) |22|I23) and later independently by Harris, Kanti 
(2005) [25] and Duffy et al. (2005) |26]. In [22l 122 , the super-radiance modes in u = u - mQ < 



are explicitly found for the first time for the MP black holes (see Fig. 14). In the winter 2005-2006, 
Casals et al. (2005) [27J found the first numerical results for vector field (s = 1) and Ida-Oda-Park 
(2006) |24] found numerical results for a vector as well as a spinor (s = 1/2, 1). Casals et al. (2006) 
[28) also got similar results for s = 1/2. Even though the effect of the mass and charge of radiated 
particles is expected to be small for semi-classical black holes, Mbh ^ M^, for which we can only 
neglect the back reactions to the black hole geometry, some related studies have been conducted 
for massive (charged) scalar emission by Casals et al. (2008) [29] . Sampaio (2009, 2009) |30| l3T] 
and Kanti et al. (2010) [32j . for bulk massive scalar by a tense brane by Kobayashi et al. (2007) 
|33) and for massive fermion by Rogatko et al. (2009) ^34j- Also a part of gravitational radiation 
(to tensor mode) was studied by Cornell et al. (2006) [35], Cardoso et al. (2006) [36j and Kanti et 
al. (2009) [37|. 



In Ida-Oda-Park (2006) [24] . the results for the time evolution of the MP black hole were found 
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FIG. 14. The power spectrum of the D = 10 dimensional black hole to the brane vector field and a 
schematic of horizon geometry. The highly rotating black hole emits Hawking radiation and settles down to 
a non-rotating black hole. When the black hole is non-rotating, the spectra are 'simple' and 'calm' compared 
to the violent and exciting radiation from a highly rotating black hole with many angular modes with several 
(I, to). By reconstruction of the spectrum, in principle, one can deduce the angular momentum and mass of 
the decaying black hole. 



and the amount of energy loss through the Hawking radiation to the complete set of the standard 
model particles (quarks, leptons, Higgs boson and gauge bosons) was calculated for the spin-down 
and Schwarzschild phases. The result is expected to be valid for a lower dimensional case (D = 5), 
especially if the Kaluza-Klein gravitons are heavy as in the Randall-Sundrum model, but a sizable 
modification is expected in higher dimensional cases {D > 6) due to a sizable gravitational radiation. 



Below is a summary of the time evolution of a rotating black hole (see Fig. 11 for a schematic of 
black hole evolution). 



Generically, the faster the black hole rotates, the larger the rate of the angular momentum 
radiation. Therefore, fast angular momentum loss happens during the early stage of evolution, 
then the phase with a rather low emission rate persists until the angular momentum is 
completely radiated away and the black hole stops rotating. During this spin-down phase, 
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typically more than 70% or 80% of a black hole's mass is lost before the angular momentum 
parameter becomes very small (< 1% of the initial value) depending on the number of extra 
dimensions. 

• The most effective loss of the angular momentum is caused by vector field radiation. Com- 
pared to the spinor and scalar field, the vector field emission dominates the radiation when a 
black hole rotates fast. Thus, the Hawking radiation is not equal for each spin state. In the 
standard model, the gluon (with eight colors) and weak gauge bosons (VF^, Z^) and photon 
(7) contribute and dominate over fermions (quarks and leptons) and the Higgs boson. 

• For a typical event with a large initial angular momentum, more than half of a black hole 
mass is emitted when the hole is still rotating such that the Hawking radiation would be 
anisotropic rather than isotropic. 

• At the end of the decay, the mass of the black hole will reach the Planck scale: M ~ 
M£). This domain is not controlled by semi-classical GR. A sizable back-reaction due to the 
discontinuous radiation of particles could also be important in this domain as the mass of 
the black hole is not large. 

A theoretical issue of separability of the Hamilton- Jacobi equation and especially the Klein- 
Gordon equation is important to predict the detailed behavior of fields in the black hole background. 
The MP black hole posseses spacetime symmetries generated by the Killing vectors and also admits 
the antisymmetric Killing- Yano and symmetric Killing tensors. The spacetime symmetry of the 
D = 2n + 1 + e dimensional MP metric guarantees the existence of n -|- 2 integrals of motion 
including the normalization of the velocity. To ensure separability, at least n = D — (n + 2) 
additional integrals of motion must exist. For D = 4 and D = 5 one has n = 1, thus the Killing 
tensor that exists in the MP metrics is sufficient for the separation of variables. For D > 6, n > 1, 
such that one cannot expect the separation of variables unless additional hidden symmetries exist 
[T8n[T82] . 

2. Small ADD black hole and RS black hole 

If a black hole is small enough (i.e., rh ^ re) the small black hole can be regarded as a higher 
dimensional black hole with 5'^"'"" event horizon topology in 4 -|- n dimensions. As long as the 
curvature radius is large compered to the size of the event horizon, there seems to be no qualitative 
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FIG. 15. A black hole on the brane recoils with bulk radiation. This process is forbidden in the presence 
of Z2 symmetry about the brane position. A black hole of mass M should be split into two smaller black 
holes with mass M/2 for each if Z2 is exact. However, the entropy of the initial black hole is 5ini = S{M) oc 
j^(2+n)/(i+n) g^^^ ^Yie entropy for the final black holes is Sfinai = 2S{M/2) ex 2(M/2)(2+»)/(i+") < Si^u 
smaller than the original entropy. Thus, such a process is forbidden by the area theorem. 



distinction between a black hole in the ADD background and the one in the RS background; both 
are well approximated by MP black solutions. However, it should be noted that RS geometry has 
a Z2 reflection symmetry about the 'brane' where the mini black hole is likely to be produced such 
that the behavior of a black hole in different scenarios might be distinguishable |159| . A small ADD 
black hole (i.e., a black hole in the ADD background) radiates mostly to the bulk in the first phase 
of production, particularly due to the large angular momentum. A black hole can recoil and leave 
the brane |167fll70] . However, for a small RS black hole (i.e., a black hole in the RS background), 
bulk radiation is strongly suppressed due to the small number of bulk degrees of freedom and the 
presence of Z2 symmetry. A black hole cannot recoil and leave the brane. 



One can easily understand such a phenomenon. Consider a black hole on a brane as in Fig. V B 2 
In the presence of Z2 in RS, recoil of the black hole is identified with splitting of the black hole into 
two smaller black holes in opposite directions. Each smaller black hole has a half the mass of the 
original black hole. Now considering the entropy of the initial and final states, one can conclude 
that the final state has a smaller entropy than the original entropy, which is forbidden by the area 
theorem (the black hole area never decreases). As the recoil and bulk radiation of graviton should be 
regarded as non-visible or missing energy signatures at particle detectors, ADD and RS black holes 
can be distinguished by the amount of missing energy signatures, for instance. However the details 
of the signature could quite possibly depend on the details of brane configuration and brane tension. 
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which should be clarified by further studies. The role of more general warped compactifications 
also remains to be understood. 

C. Recent CMS results 

Searches for microscopic black hole production and decay in pp collisions at a CM energy 7 TeV 
have been conducted and reported by the LHC (CMS collaboration) using 35 pb~ |161| and 1 fb^ 
|162| (CMS physics analysis summary (CMS PAS EXO-11-071) ) data sets obtained before the 
summer of 2011. More data has been collected by the CMS as well as ATLAS so that new results 
are expected soon. The higher energy upgrade from 7 TeV to 14 TeV of the LHC is expected, and 
we will definitely will be able to learn more about the physics on TeV scales in the coming years. 
The higher energy upgrade is particularly important for the discovery of the new physics models 
with a high energy threshold. 

For the black hole search at the LHC, the following set of criteria and assumptions on signal 
and background has been considered in the CMS experiment |16H I162J : 

• The ADD model is assumed for low scale gravity with Md ~ 1 TeV. The produced black 
hole is confined on a brane where all particles except gravitons reside. 

• A semi-classical black hole has been considered even though a large quantum correction is 
expected (we will discuss this more in the later part of this section) . 

• High multiplicity (A^ > 2 — 8) from the black hole with a large entropy and democratic, 
highly isotropic Hawking radiation, which is typical of a non-rotating black hole, is taken 
into account for signals. 

• For a trigger, the final state particles (mostly jets and hadrons) carrying more than 100 GeV 
of energy are considered. The large scalar sum of the transverse energy, Ht = Y^jet-^T > 
50(100-200) GeV, with jet Et > 10(20-30) GeV at the first level trigger (High level trigger) 
is used. For event selection, a large scalar sum of the Et on the A^ individual objects including 
jets, electrons, photons and muons, (= ^^ >50GeV "^T" ^ (1-5 — 4.0) TeV), are considered 
for various Md G (1.5,3.0) TeV in D G [6,10]. (note: Md = Md x (87r/A/'„)^/("+2) and 
A/; = (27r)" in the PDG unit.) 

• Black hole signal events are simulated using the parton level BlackMax generator [38] (version 
2.01.03), followed by a parton showering fragmentation with PYTHIA |163] (version 6.420) and 
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a well validated simulation of the CMS detector response. 

The main background to black hole signals are from QCD multi-jet events, W jZ + jets 
and tt production were estimated from Monte-Carlo simulations using MADGRAPH |164) with 
the CTEQ6L PDF set |133| followed by PYTHIA parton showering and full CMS detector 
simulation via GEANT4 [165j . 



The main results are depicted in Fig. 16 and Fig. 17 (both are from |161) ). Similar but slightly 



improved results with a larger data sample are found in |162) . In summary, no significant deviation 
from the background was found in the 35 fb^ nor in the 1 fb~ data samples. Only a limit on black 



hole mass (see Fig. 17) for given values of Md in the TeV range could be set. In [162j . results for 



higher multiplicities up to A^ > 8 are presented, but still no significant deviation has been found. 
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FIG. 16. Total transverse energy St for events with multiplicites 71 > 3, n > 4 and n > 5, respec- 
tively, obtained by the CMS using data sets with 7 TeV and 35fb^^. The shaded band is the background 
(QCD+leptons+tf) for prediction with its uncertainty. The figures are from |161| . Similar results can be 
found in [TB51 for Ifb"^ data. 



Some comments on the recent CMS results are in order. As was emphasized previously |16HI162] . 
the results should be treated as indicative, rather than precise, as the semi-classical approximation 
is not valid in the parameter region covered by the experiment so far |166) . There are two main 
cautions in interpreting the results: 



The parameter range in the CMS results is actually out of the validity range of the semi- 
classical approximation of a black hole. 



Interpretation of the results depends on particular realization of the low energy gravity mod- 
els. 
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FIG. 17. The 95% confidence level upper limits on the produced black hole cross-section and mass in the 
CMS results. The figures are from |161| . Similar results can be found in |162| . 



Having a limited collision energy y^ = 7 TeV, the considered 'black hole' cannot be heavy 
enough beyond the quantum gravity regime in 0(1) TeV. The ratio of the minimum black hole 
mass and the fundamental scale, which was excluded by the CMS data, is in the range |16H 1166] 
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FIG. 18. The higher order curvature term, A, (top) and entropy (bottom) are plotted for a higher 
dimensional black hole {D — 6(dotted), D = lO(solid)). The vertical column in gray is the CMS exclusion 
region: M/Md = [3.5/3.0(4/3.5), 4.0/1.5(4.5/1.5)] for D = 6(10). 
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As was discussed in earlier sections, the validity of semi-classical approximation is guaranteed 
when the black hole is massive enough, M ^ M/j. However, in the claimed exclusion range, the 
large quantum gravitational correction of the order of A = ( — 7^^ J > 100% with a positive 
(unknown) power p > is expected. We may not impose any geometrical meaning on an object 
with a small entropy. See an estimation of higher order correction (upper) and the extrapolated 



value of entropy (lower) in Fig. 18 for D = Q (dotted lines) and D = 10 (solid lines). For consistency 
of the perturbative expansion, the correction term must be significantly smaller than the leading 
order term (A ^ 1) and entropy should be large {S ^ 1), but in the CMS exclusion range (vertical 
columns in gray scales), the correction is expected to be as large as the leading term. The entropy 
of the black hole is 0{\). Within this parameter space, all the MC simulations suffer from large 
quantum corrections, which can lead a significant change in the final result. There may be a chance 
for the LHC to form some black hole precursors |127| I177J and low energy Kaluza-Klein excitation 
modes of bulk fields]^ However, a semi-classical black hole now seems out of the reach of the LHC 
run with ^/s = 7 TeV. 

D. A guide for future improvement 

Here is a guide for future improvement |171) . 

• The detailed brane realization of the standard model is still unknown (see Sec. [H] ). 

• The balding phase should be understood through dynamical simulation, which can most 



probably only be done via a numerical study (see Sec. IVB3|. 



For D ^ 4, spin-2 graviton emission to the bulk from highly rotating black holes can be 
sizable and dominant. Hawking radiation to a spin-2 particle has been considered only for 
the non-rotating case in higher dimensions D > 4 |35fl37) . 



The final state of a black hole, i.e., the Planck phase, is extremely poorly understood. Full 
quantum gravitational consideration is required (see e.g. |172fll76] for some recent discus- 
sions). 

Details of signals should depend on many factors (e.g., 2 to 2 dominance |177| . inelastic 
effect [178], recoil [THTHTTn] . split-brane [I79l [ED], hidden symmetry [ml [l82], instabihty 
in relation with quasi- normal mode and black hole bomb |183l 1184) and |185t4187] . etc.). 



An important counter example is the model in the presence of a compact hyperbolic extra dimension background 



where no low energy excited Kaluza-Klein modes are expected (see Sec. II D I 
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• The quantum gravity correction to black hole geometry and Hawking radiation should be 
modeled; see |188) (Gauss-Bonnet corrected Black hole) and |189j (non-commutative geome- 
try inspired black hole). 

VI. CONCLUSIONS 

The LHC is the most promising place to look for new physics models on the TeV scale including 
the low scale gravity models such as ADD and RS. If one of the low scale gravity models is realized 
in nature, the microscopic black holes on the 1 TeV scale will be produced at the LHC and their 
signatures, especially the Hawking radiation, might provide a valuable guide to providing a deeper 
level of understanding of quantum gravity and spacetime structure. In this article, we review the 
progress of understanding this exciting area of research over the past decade, emphasizing both 
theoretical and experimental efforts. On the theoretical side, analytic as well as numerical studies 
have been performed to understand the dynamical process of black hole formation in high energy 
collisions mostly in the semi-classical domain. The Hawking radiation to brane localized fields (and 
also partly to the bulk fields) from a higher dimensional rotating black hole is understood via the 
precise calculation of greybody factors. Monte-Carlo simulation codes have been developed for the 
"black hole events" at the LHC. On the experimental side, the LHC begins to produce data for the 
TeV scale gravity models and will reach the true trans-Planckian or microscopic black hole domain 
in the coming years after its collision energy is upgraded to 14 TeV from the current 7 TeV as the 
LHC was originally designed. Even though the current understanding is limited due to the lack of 
a full quantum gravity model, we hope to learn more by observing real data from the TeV domain. 

Finally, a short comment on the safety of the LHC is in order. The recent analysis shows that the 
macroscopic effects of TeV-scale black holes should have already been seen in various astrophysical 
environments if those black holes are dangerously long-lived or even stable |14| [TS] . This analysis 
is important since it excludes the possibility that any mini black holes that might be produced by 
the LHC could cause catastrophic damage. If the LHC can produce black holes, there are various 
places in the universe such as neutron stars and white dwarfs where similar mini black holes can 
also be produced at an even more frequent rate by collisions of the ultra-high energy cosmic ray 
particles and nucleons in those dense astronomical objects. Basically all of those dense astrophysical 
objects are of the same or a higher level of danger. Thus, the null observed signals from those dense 
astrophysical objects ensures that there is no risk of any significance whatsoever from such black 
holes. 
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Appendix A: Appendix: Planck units in D dimensions 



The Planck units in Z) = 4 + n dimensions can be determined as: 
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[h] = l'^t-^m, 

[c]=LT-\ 



G%h'^c' 



'D 



1 



fGDh\T^ 



Md 



Gd 



(A.l) 
(A.2) 
(A.3) 
(A.4) 

(A.5) 
(A.6) 
(A.7) 



The Schwarzschild-Tangherhni black radius (r^) corresponding to the mass M and impact pa- 
rameter of the scattering b in terms of angular momentum: 



Ts = Cd 



/GdM\D-3 



6 = 2 



J 

7^' 



(A.8) 
(A.9) 



where Cd = ^0-^0-2 ) " ^ ^^^ ^^^2 = raD-syl) - 



At the limit M = ^ > Md with h^O and fixed Gd, 



iz 


) — > 


0, 


AdB = 


Vs 


Ts 


— > 


finite, 



(A.IO) 
(A.ll) 
(A.12) 



such that one can conclude that in the trans-Planckian domain, the dynamics would be described 
by rs >^D,AdB- 
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Appendix B: Appendix: Solving the generalized Teukolsky equation 



The solution to the generahzed Teukolsky equation can be found numerically |21H24) . To get 
the solution, one should impose the 'purely ingoing wave' condition near the horizon (r — t- r^) then 
numerically integrate the differential equation to the far field regime (r — )• oo) where ingoing and 



outgoing waves are superposed (see Fig. 19). 



Near Horizon 



Far Field 




Numerical integration from NH to FF 



NH 
■ingoing 

purely ingoing wave 



ingoing + outgoing waves 



FIG. 19. Schematic summary of the method to find the greybody factor by solving the generalized Teukolsky 
equation (details in |24|V 



On the (3 + 1) dimensional 'brane', the induced metric of the (4 + n)-dimensional Myers-Perry 
solution with a single nonzero angular momentum (on the brane) is given by 



ds^ 



A - a^ sin^ •& , , 2a(r'^ + a^ - A) sin^ i} , , 
dt^ + -^ '- dtd^ 

S ^ A 



(B.l) 



where 



S = r^ + o^ cos^ ??, A = r^ + a^ - ^r^-". 



(B.2) 



The parameters fi and a are equivalent to the total mass M and the angular momentum J evaluated 



at the spatial infinity of the (4 + n)-dimensional space-time (see Sec. III). 

The radial equation of the generalized Teukolsky equation is given by a linear, second-order 
equation as shown in |21] : 

d^R 



ar-^ dr 



(B.3) 
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where 7] and r are determined in the Kerr-Newman frame as: 

(s-l)A' + 2iK 2iu;r(2s - 1) - \ ,„ ^, 

ri = ^ ,r= (B.4) 

where K = (r^ + a'^)w — ma. 

1. Near horizon: r -^ r^ 

To find the near horizon (NH) solution, it is convenient to expand all quantities in the (r — Tj^)^ 
series: 



A = ^Afc(r-r;,)^ (B.5) 

fc=i 

K = Y,Kk{r-rh)\ (B.6) 

fc=i 

^=Y^ rik{r-rhf, (B.7) 

k=-l 

r=Yl ^k{r-rhf. (B.8) 



fc=-i 
Then the solutions near the horizon (NH) r ^ r^ are read 

<« = l+ai(r-r,) + |(r-r,)2 + ..., 

<'? = e'*"'*(r-r,)'^(l + 6i(r-r,) + •••), (B.9) 

where the expansion coefficients Oj's and 6j's are found as follows: 

ai = - — , (B.IO) 

(r?o + r_i)ai+ro 
a2 = ^^^^ , (B.ll) 

fB.12) 



The problem now is to integrate Eq. (B.3) from purely ingoing initial conditions at r = r^ out 
to r — 7- oo. However, i?ou{ is unstable against contamination of the purely ingoing solution. To 
counteract the above contamination, let 



R = R-{1 + ai{r - r^)) « ^(r - rnf + ■■■ (B.13) 
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in the vicinity of the horizon, which satisfies the following equation: 

CR = g, (B.14) 

where C = cP / dr'^ + rjd/ dr + t and g = — £(l + ai(r — r//)) = —rjai — t (1 + ai(r — r//)). Equation 



(B.14) can be integrated numerically from the near horizon to the far field regime. Calculation 



details could be found in |[24j. 



2. Far field limit : r ^ oo 

At the far field limit (FF), the solutions to the radial equation can be found easily. 
For a spinor field with spin s = 1/2, 



Iff 

fi-i ^•1 /^^ 



-]f jf 

' 9 



<u"t~e2^'=^4(l + - + %•••), (B.15) 



where 



(BAT) 



Thus i?i /2 is read 



Ry2{r ^ oo) ~ (Fin + 1 - ai) + air + -^ + l^out^^- (B.18) 



Using this expression, we can easily read out Ys without numerical difficulties at FF. The 
greybody factor is found as follows: 

9, > W |2 

rs=i/2 = i-^^. (B.19) 



5=1/2 - ^- 777 1^ |2 • 



For a vector field with spin s = 1, 



Co Co 






rt* ^ /y* fy*^ /y^'j 



where 



2uj 
A^ — Aau:{aijj — m) 



Then, Ri is read 



Jliiur^ 
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Rll-e^'''^\l+'^l + % + % + ---l (B.20) 



c\ = -i^, (B.21) 



C2 ^-2 ' ^^-22) 



R^(^r ^ oo) ~ (l^n - ai) r + (YinC^' - 1 + ai) + ^^ + Fout ^^. (B.23) 

r r 

Similarly when s = 1/2, one can easily read out Ys without numerical difficulties at FF. Finally, 
the greybody factor for a vector field can be determined as 

9/ i2 \V |2 

r.-i = l-FV^^. (B.24) 

L,D v: 2 ^ ■' 

I <-2 1 I 111 I 
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